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We present the first upper-limit measurement of the branching fractions of the isospin-violating 
transitions Xc 0,2 —>■ 'x^'Oc. The measurements are performed using 106 x 10® -!/>(3686) events accumu¬ 
lated with the BESIII detector at the BEPCII e''"e“ collider at a center-of-mass energy corresponding 
to the if (3686) mass. We obtained upper limits on the branching fractions at a 90% confidence level 
of B{xcO —> Ti^Vc) < 1.6 X 10“® and B{xc2 —>■ 7r°»?c) < 3.2 x 10“®. 

PACS numbers: 13.25.Gv, 13.20.Gd. 


I. INTRODUCTION 


rates of isospin-symmetry breaking modes are in gen- 


Isospin is known to be a good symmetry in the 
hadronic decays of charmonium states. The decay 
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eral found to be very small. For example, the branch¬ 
ing fraction {B) of the measured isospin-violating tran¬ 
sition tf}{3686) —?> TT^J/ip was found to be only (1.26 ± 
0.02{stat.) ± 0.03{syst.)) x 10“^ [T], whereas for other 
hadronic transitions such as •0(3686) —>■ 7r+7r“J/0, the 
branching fraction is (34.45 ± 0.30) x 10“^ [2] and thus 
significantly stronger. 

Although isospin breaking is found to be very small 
for the conventional charmonium states, the mysteri¬ 
ous X(3872) resonance above the DD threshold decays 
strongly via the transition A'(3872) —>■ tt+tt” J/0, where 
the invariant-mass spectrum of the tt+tt” pair shows a 
clear p signature HE] and, hence, is compatible with an 
isospin-violating decay. A possible interpretation is that 
the Ar(3872) is a molecular state composed of a bound 
D*^-D^ meson pair imm)- Such an explanation is par¬ 
ticularly popular, since the mass of the X(3872) is close 
to the sum of the and masses, pointing to a 
state that could be weakly bound by the exchange of 
a color-neutral meson, similar to the deuteron. More¬ 
over, in such a scenario, the strong isospin-breaking de¬ 
cay rate of the X(3872) might be explained by the large 
mass gap between the D*^ — and the D*~^ — D~ 
(Z1+ — D*~) thresholds [TT] . A better understanding 
of the isospin-breaking mechanism in a complementary 
and well-established charmonium system below the open- 
charm threshold could be crucial to shed light on the 
nature of the X(3872). 

On the quark level, the isospin-symmetry is broken 
due to the electromagnetic interaction and due to differ¬ 
ences in the up- and down-quark masses (to„ and mo)- 
It is, therefore, believed that isospin-breaking decays can 
be used to access the up- and down-quark mass differ¬ 
ences once the electromagnetic effect is either well un¬ 
derstood or found to be negligible. An example observ¬ 
able that has been proposed to obtain the quark mass 
ratio, mu/rud, is a measurement of the ratio between the 
branching fractions of the transitions •0(3686) —)■ 7r° J/0 
and ^(3686) —)■ ?7J/0. Based on a leading-order QCD 
multipole expansion |12j and the BESIII measurement 
of this ratio [T], the up-down quark mass ratio is ex¬ 
tracted to be mu/md = 0.407 ± 0.006. This result is 
smaller than the result, rriu/rnd = 0.56, obtained using 
the Goldstone boson masses from a leading-order chiral- 
perturbation theory |13j . It is important to understand 
such a large discrepancy between the values of m„/md 
obtained on the basis of different theoretical conjectures. 

The most promising developments in this field are 
based upon an effective-field theoretical approach. A 
nonrelativistic effective-field theoretical (NREFT) study 
by the Jiilich and IHEP groups suggests that inter¬ 
mediate (virtual) charmed-meson loops are the domi¬ 
nant source for the isospin breaking in the transition 
0(3686) —•7r°J/0 [T3J[1^. According to the proposed 
theory, the contribution of charmed-meson loops to the 
amplitude of the process is enhanced by a factor of 
{v/c)~^ ~ 2, where v is the heavy-meson velocity in the 
loops. Detailed studies of different isospin-violating tran¬ 


sitions in charmonium below the DD threshold and the 
effect of virtual charmed-meson loops on the widths of 
the transitions are described in Ref. [TB] . 

The NREFT calculations described above are based on 
a first estimate, exploiting diagrams involving the lowest- 
lying pseudoscalar and vector charmed mesons following 
heavy-quark symmetry and chiral symmetry. Although 
these theoretical calculations give qualitative insights in 
the isospin-breaking mechanisms in charmonium decays, 
the authors in Ref. m state that only with a further 
developed effective-field theory that includes Goldstone 
bosons, charmonia, and charmed mesons as the degrees 
of freedom, it would be possible in the future to ex¬ 
tract the light-quark masses from quarkonia decays. Cur¬ 
rently, for such a theory, quantitative predictions of in¬ 
dividual branching fractions of isospin-forbidden decays 
of charmonium are difficult, because information on the 
coupling constants fpoD between different charmonium 
states and DD-mesons is limited. The theory requires 
constraints from experimental data, in particular from 
measurements of decay rates of other isospin-violating 
transitions in charmonium m- 

In this paper, we present an experimental study of the 
isospin-suppressed transition of the charmonium P-wave 
states Xc 0,2 to the ground state rjc via the emission of 
the TT®. The Xc 0,2 states are obtained via electromagnetic 
transitions, 0(3686) —> 7 Xc 0 . 2 , whereby the 0(3686) reso¬ 
nance is directly populated via the e+e“ annihilation pro¬ 
cess. The transition Xci is not considered in this 

analysis since it violates conservation of parity and angu¬ 
lar momentum. According to Ref. |16] , the dimensionless 
suppression factor for the loops in Xco —t is 0.2. This 

factor is smaller than in the process 0(3686) —>■ 7r°J/0, 
however, through the interference with the tree-level am¬ 
plitude, meson loops may give a significant contribution 
and cannot be neglected. 


II. THE BESIII EXPERIMENT AND DATA SET 

The analysis is based on the 0(3686) data sample ac¬ 
cumulated by the BESIII detector in 2009. The total 
number of 0(3686) events is (106.41 ± 0.86) x 10® |17j . 
corresponding to an integrated luminosity of 156.4 pb“^. 
In addition, 42.6 pb“^ data collected at a center-of-mass 
energy of 3.65 GeV, are used to estimate the background 
from nonresonant processes. 

The BEijing Spectrometer III (BESIII), described in 
detail in Ref. m, is a detector for r-charm studies run¬ 
ning at the Beijing Electron-Positron Collider (BEPCII). 
BEPCII is a double-ring e+e” collider with a designed 
peak luminosity of 10®® cm“^s“® at a beam current of 
0.93 A. The cylindrical core of the BESIII detector con¬ 
sists of a main drift chamber (MDC), a plastic scintil¬ 
lator time-of-flight system (TOF), and an electromag¬ 
netic calorimeter (EMC), which are enclosed in a super¬ 
conducting solenoidal magnet providing a 1 T magnetic 
field. The solenoid is supported by an octagonal flux- 
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return yoke with resistive-plate chambers forming a muon 
counter system. The MDC is a small-cell, helium-based 
(40% He, 60% CaHg) subdetector consisting of 43 layers 
and providing an average single-hit resolution of 135 /rm, 
and a charged-particle momentum resolution of 0.5% at 
1 GeV/c. The EMC subdetector consists of 6240 CsI(Tl) 
crystals in a cylindrical structure (barrel) and two end 
caps. For 1 GeV photons, the energy resolution is 2.5% 
(5%) and the position resolution is 6 mm (9 mm) for 
the barrel (end caps). The TOE system consists of 5 cm 
thick scintillators, with 176 detectors of 2.4 m length in 
two layers in the barrel and 96 fan-shaped detectors in 
the end caps. The barrel (end cap) time resolution of 
80 ps (110 ps) provides 2a K/tt separation for momenta 
up to 1 GeV. 

To optimize the event selection, to estimate back¬ 
ground contributions, and to evaluate the detection ef¬ 
ficiencies, Monte Carlo (MG) simulated samples are 
obtained exploiting a realistic model of the detec¬ 
tor. For this, the GEANT4-based simulation software 
BOOST [19] is used which includes the geometry and 
material description of the BESIII spectrometer, and the 
detector response. A MG sample based on 106 M in¬ 
clusive '0(3686) decays is used to study the background. 
This inclusive sample is generated with KKMC |20| plus 
EvtGen [211 [22] and the known branching ratios are taken 
from the Particle Data Group (PDG) [5], while the un¬ 
known ratios are generated according to the Lundcharm 
model [23]. The decay modes Xc 0,2 —t are not 

present in the inclusive MC simulation. Signal MG sam¬ 
ples are generated to determine the detection efficiency 
and to model the signal shape. In the MC simulations for 
the processes presented here, the '0(3686) —>■ xXcJ decay 
is assumed to be a pure El transition, and the polar an¬ 
gle, 0, follows a distribution of the form 1 -|- a cos^ 9, with 
0 = 1 and 1/13 for J = 0 and 2, respectively [Ml 1^ . 
The Xc 0,2 and rjc —>■ KgK^n^ decays are as¬ 

sumed to be pure phase-space decays. 


III. DATA ANALYSIS 

For the identification and selection of '0(3686) —>■ 
7XcO,2 —t events, where 7r° — >■ 77, t]c — >■ 

the Kg is reconstructed in its decay mode to 7r“'"7r“, re¬ 
sulting in the final state SydTrA, where K and tt are 
charged. 


A. Event selection 

Charged tracks are reconstructed from the MDC hits. 
For each charged-particle track, its polar angle must 
satisfy |cos0| < 0.93. A good charged-particle track 
(excluding those coming from a Kg) is required to be 
within 1 cm of the e“''e“ annihilation interaction point 
(IP), transverse to the beam line and within 10 cm of 
the IP along the beam axis. Charged-particle identi¬ 


fication (PID) is based on combining the energy loss, 
dE/dx, in the MDC and TOF information to construct 
PID chi-squared values X^ioi'^)^ ^hat are calculated for 
each charged-particle track for each particle hypothesis i 
(pion, kaon). 

Photons are reconstructed from isolated showers in the 
EMC. The showers in the angular range between the bar¬ 
rel (|cos0| < 0.8) and end caps (0.86 < |cos0| < 0.92) 
are poorly reconstructed and excluded from the analysis. 
Good photon candidates must have a minimum energy 
of 25 (50) MeV in the barrel (end cap) regions. EMC 
timing requirements are used to further suppress noise 
and energy depositions unrelated to the event. 

Events with four charged-particle tracks with a net 
charge of zero and at least three good photon candidates 
are retained for further analysis. 

Kg candidates are reconstructed from secondary ver¬ 
tex fits to all the charged-track pairs in an event (with 
a pion-mass assumption). Candidates with an invari¬ 
ant mass within 10 MeV/c^ of the Kg nominal mass 
are considered and the combination with the smallest 
chi-squared of the vertex fit is chosen. The event is 
kept for further analysis if the secondary vertex is at 
least 0.5 cm away from the IP. The reconstructed four- 
momenta of the 7r“'" and 7r“, corresponding to the Kg 
decay, are used as input for the subsequent kinematic 
fit. To suppress the KgKg background, the remaining 
charged-particle tracks are required to not form a good 
Kg candidate. The 7r° candidates are reconstructed from 
pairs of photons with the invariant mass in the range 
0.11 < Mjy/{GeV/c^) < 0.16, with the resolution 
of about 5 MeV/c^. 

The S^SttK candidates are then subjected to a four- 
constraint (4C) kinematic fit, with the constraints 
provided by four-momentum conservation. The dis¬ 
crimination of charge-conjugate channels {KgK~^Tr~ or 
K^K- '7r+) and the selection of the best photon candidate 
of the 0(3686) —>■ xXcO ,2 transition among multiple candi¬ 
dates are achieved by taking the event with the minimum 
= xlc + Xhoi^) + XpidM’ where xlc is the chi- 
squared of the 4C kinematic fit. The 7r° is reconstructed 
from the two-photon combination with an invariant mass 
closest to that of a neutral pion. Events with xic < 
and with an invariant mass of the reconstructed rfc, 
^KOK± 7 rT, in the range 2.70 < /(GeV/c^) < 

3.30 are accepted for further analysis. The maximum 
value of X 4 C is determined by optimizing the statisti¬ 
cal significance S/\/S + B in the rjc signal region, where 
S (B) is the number of signal (background) events ob¬ 
tained from the signal (inclusive) MC samples. For the 
estimate of S, the branching fractions of Xc 0,2 —t are 

assumed to be 10“^ in analogy with the isospin-violating 
process 0(3686) —>■ [I]- The signal region is de¬ 

fined as 2.90 < M^o^±^^/(GeV/c^) < 3.05. The Xco 
and Xc 2 signal regions are defined for transition-photon 
candidates with energies in the KgK^n^ center-of- 
mass system, E^, in the ranges of 0.24 < E^/{GeXf) < 
0.28 and 0.10 < E^/{GeXf) < 0.15, respectively. 
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FIG. 1. Invariant-mass distributions of candidates (a and b) and K^tv^ candidates (c and d) for XcO (a and c) and Xc 2 

(b and d) mass regions, respectively. Dots represent data, filled histograms represent inclusive MC results, open histograms 
show results of signal MC simulations based on a phase-space distribution (arbitrarily scaled). 



FIG. 2. Invariant-mass distributions of . Left: XcO —>■ right: Xc 2 —^ Vc . Dots represent data, 

filled histograms represent the inclusive (light) and arbitrarily scaled signal (dark) MC results. The peak around 3.12 GeV/c^ 
is due to the ip{3686) —>■ Tv^-jr^J/ip background channel. 


FIG.[^shows the invariant-mass distributions of 
candidates [(a) and (b)] and candidates [(c) and 

(d)] with ir'^KgK^TT^ masses within the XcO [(a) and 
(c)] and the Xc 2 [(b) and (d)] mass regions. The most 
prominent peak (with the highest intensity and narrow¬ 
est width) stems from decays involving a K*{892). Those 
are evidently background processes, because the chan¬ 
nel of interest, Xc 0,2 —^ with t 

cannot involve K*{892)^^^^ —>■ decays, since 

the latter does not involve a 7r°. The regions 0.84 < 
M^o(±)^o/{GeV/c^) < 0.95 are excluded in the further 

analysis to suppress the background from iir*(892)° and 
iL*(892)^ decays. This condition is optimized to obtain 
the best statistical significance of the signal. 

FIG. [2] shows the invariant-mass distributions of 
KgKir for candidate events with masses corre¬ 
sponding to XcO ,2 —t transitions. The data show no 

visible peak in the rjc signal region. In the present anal¬ 
ysis, upper limits at the 90% confidence level (C.L.) for 


the transitions Xco.2 are determined. Inclusive 

MG results do not reproduce the number of events found 
in the data, but reproduce the shape of the invariant- 
mass distributions of K^Kir quite well. The discrepan¬ 
cies between data and inclusive MG are primarily due to 
inaccuracies of branching fractions in the generator and 
due to mismatches in their corresponding decay dynam¬ 
ics. We note, however, that the inclusive MC data have 
solely been used to identify background sources and to 
optimize selection criteria. The selection criteria were op¬ 
timized by assuming a branching fraction of 10“^ for the 
Xcj channels in combination with background 

taken from the inclusive MC sample. We varied the 
signal-to-background ratio by a factor of 2, adjusted our 
selection criteria accordingly, and found a negligible ef¬ 
fect on the precision of our final result. 

Efficiencies are calculated using the signal MC simula¬ 
tion samples and are found to be 5.8% and 8.6% for the 
XcO and Xc 2 channels, respectively. 
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B. Background studies 

Background events from ^/^(SGSd) decays are studied 
with the inclusive MC sample. These studies showed 
that the channel '!/'(3686) —?> TT^n^J/ip, J/i/^ —>■ KgK^ir^ 
results in a peak around 3.12 GeV/c^ in the KgKir 
invariant-mass spectrum as can be observed from Fig. 

In this type of transition, one of the photons originating 
from TT^ decays may escape, which causes a smaller total 
energy for the event. The kinematic fit increases the en¬ 
ergy of the charged-particle tracks, which results in a shift 
in the invariant mass from 3.10 GeV/c^ to 3.12 GeV/c^. 
This decay channel is taken into account in the final fit 
to the invariant-mass spectrum of KgKir as described 
below. 

The major background contribution stems from the 
channels i/;(3686) —>■ "jXcJ, XcJ —^ n^KgK^Tr^. These 
channels have final states that are kinematically iden¬ 
tical to the signal of interest and, therefore, cannot be 
removed easily. Partly, this type of background has been 
suppressed by vetoing iF*(892)° signals via a cut on 
the Kg’^TT^ mass since the background decay, Xc 0,2 —> 
ir'^KgK^TT^, contains intermediate K*{892) resonances, 
as discussed earlier. We note that the remaining contri¬ 
bution of this type does not result in a peaking back¬ 
ground in the signal region. 

The background contribution from e“'"e“ —^ // pro¬ 
cesses, where / = e,/r, d, u, s, is studied using the con¬ 
tinuum data taken at ^/s = 3.65 GeV, and it is found 
to be negligible. Using the exclusive MG simulations 
and taking the corresponding branching fractions from 
the PDG [2], the contribution of XcJ n^KgK^Tr^ 
and tfj{3686) —)■ 7r°7r°J/'0 channels in the region 2.70 < 
MKOKTr/{GeV/c^) < 3.30 is found to be 2260±340 and 
1668±260 events for the Xco and Xc 2 selection criteria, 
respectively, where the errors are mainly due to the un¬ 
certainties in the branching fractions. The total num¬ 
ber of data events in the same region is 2477±50 and 
1527±39, respectively. These are compatible within the 
uncertainties. No significant peaks are observed in the 
signal region. 

C. Upper limits for the number of signal events 

To extract the number of rjc events, an unbinned max¬ 
imum likelihood fit is applied to the candidate events 
with KgK^TT^ invariant-mass distributions in the region 
2.70 < Mj^oj^^/{GeV/c^) < 3.30. The rjc signal is de¬ 
scribed by a Voigtian function, which is a Breit-Wigner 
function convoluted with the detector resolution. Param¬ 
eters of the Breit-Wigner function are taken from the 
PDG [2], and the detector resolution is obtained from 
a fit to the signal MC set. These parameters are fixed 
while fitting the data. From the background studies, no 
peaking background is expected in the signal region. The 
smooth background is described by a third-order Cheby- 
shev polynomial. A Voigtian function and a Landau 


plus Gaussian function are used to describe the structure 
around 3.12 GeV/c^ for the XcO and Xc 2 mass regions, 
respectively. The line-shape parameters of the structure 
around 3.12 GeV/c^ (for both Voigtian and Landau -I- 
Gaussian functions) are fixed to the values obtained from 
the exclusive MC sample. The MC sample was obtained 
by simulating the channel ip{3686) —>■ 7r°7r°J/'0 with the 
exclusive decay J/ijj —>■ KgRir. The total fit results are 
shown in FIG.|^ Using the maximum likelihood method, 
the upper limits on the number of signal events, , at 
the 90% C.L. are found to be 14.1 and 35.9 events for the 
XcO and Xc 2 mass regions, respectively. 


IV. SYSTEMATIC UNCERTAINTIES 

Table |l] summarizes all the systematic uncertainties 
that are considered in the analysis. Below we discuss 
in more detail the individual sources and the procedure 
that is used to estimate the errors. 

The tracking efficiency for kaons as a function of trans¬ 
verse momentum has been studied using the process 
J/ij} —>■ , Kg —>■ 7r+7r“ and the tracking effi¬ 

ciency for pions (not originating from Kg) as a function 
of transverse momentum has been studied using the pro¬ 
cess ■!/’(3686) —>■ 7r+7r“ The difference in efficiencies 
between data and MC simulations is 2% for each K or 
TT track. This value is taken as the uncertainty in the 
tracking efficiency. The systematic uncertainty due to 
the Kg reconstruction is 4.0%, as reported in Ref. [SH] . 
The uncertainty in the photon reconstruction is taken as 
1% per photon as reported in Ref. m- In this analysis, 
there are in total three photons in the final state, which 
yields a total systematic uncertainty due to the photon 
reconstruction of 3%. 

Some differences are observed for the xlc distribu¬ 
tions between data and MC simulations. These differ¬ 
ences are mainly due to inconsistencies in the charged- 
track parameters between data and MC simulations. We 
apply correction factors for various K (tt) track param¬ 
eters that are obtained from the control data samples 
J/ijj —> (f)TT'^TT~, (p —)■ K'^K~. The correction factors 
are used for smearing the MC simulation output, so that 
the pull distributions properly describe those of the ex¬ 
perimental data. Differences between the detection effi¬ 
ciencies obtained using MC simulations with and without 
these corrections are taken as an estimate for the corre¬ 
sponding systematic uncertainties. The uncertainties are 
1.2% and 0.8% for the XcO and Xc 2 selection conditions, 
respectively. 

A phase-space (PHSP) model used for MC genera¬ 
tion of 77c —>■ KgK^TT^ events does not include possi¬ 
ble intermediate resonances between the final-state par¬ 
ticles, for example, 1 X 1 ( 1430 )°’^. These resonances 
are observed in the test sample ip{3686) —>■ 'yr]c,r]c —>■ 
KgK^TT^. To account for K* resonances, additional 
MC samples of 7/>(3686) ^ 7Xco.2,Xco,2 ^ ^ 

Ar|(1430)iX, Ar|(1430) —>■ Kn are generated and the re- 
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FIG. 3. Fit to the invariant-mass distributions of . Left: XcO —>■ right: Xc2 —>■ 7r°?7c; Vc KgK'^'iv^. Dots 

represent data, lines represent the fit. Arrows indicate the peaks of interest. 


construction efficiencies are calculated. The difference 
between efficiencies obtained with two different genera¬ 
tor models is taken as a systematic uncertainty. For the 
Aco and Xc 2 —t decays, the corresponding 

systematic uncertainties are 1.6% and 2.7%, respectively. 

The selection of exactly four charged-particle tracks 
before the vertex cuts can introduce an additional sys¬ 
tematic uncertainty in the efficiency determination due to 
the presence of fake tracks from misreconstruction. This 
uncertainty is estimated using a sample of ^(3686) —>■ 
7 r+ 7 r“ J/i/j, J/ip —>■ 777 c, r]c decays, where 

events with at least six charged-particle tracks are ac¬ 
cepted. The fraction of events with more than six 
tracks compared to the number of events with exactly 
six charged-particle tracks are obtained for the data and 
for a corresponding signal MC sample, and the difference 
in the fractions is found to be 1%. This we take as the 
systematic uncertainty due to the preselection of exactly 
four charged-particle tracks. 

To estimate the systematic uncertainties due to the 
resolution of the transition photon, a smearing of the 
energy resolution by 1.5 MeV and a shift of 0.5 MeV 
are introduced on the MC data. The smearing gives a 
minimum when comparing data and MC line shapes, 
where the dominant contribution to the data originates 
from the processes '0(3686) — IXcJ^ XcJ —t TT^KgK^n^. 
The corresponding detection efficiencies of the channels 
of interest with the standard selection criteria have been 
calculated. The largest difference between the efficiencies 
with and without smearing for the Xc 0,2 is 0 . 6 %, which we 
quote as the systematic uncertainty due to the transition- 
photon resolution. 

The systematic uncertainty due to fitting consists of 
four parts: uncertainties due to the fitting range, the 
J/0-related background shape, the Xc 0,2 —>■ n^KgKn 
background shape, and the signal shape. The upper lim¬ 
its for the 2.70 < Mj^oj^^KGeV/c^) < 3.30 fitting range 
with a third-order Chebyshev function and fixed J/0 pa¬ 
rameters are taken as the nominal upper limit, . By 
varying the fitting ranges (2.60 < M^o^,r/(GeV/c^) < 
3.30 and 2.70 < M^oj^j^/{GeV/c^) < 3.40 for the Xco 
and Xc 2 rnass regions, respectively), we obtained a set 
of upper limits from which we take the maximum dif¬ 


ference as a systematic uncertainty. The uncertainty is 
found to be 5.8% and 15.0% for the XcO and Xc 2 mass re¬ 
gions, respectively. By changing the parameters for the 
J/'0 background from fixed to free, but using the nominal 
fitting range and the nominal order of Chebyshev poly¬ 
nomial, the result of the fitting procedure is obtained, 
and the relative difference with the nominal upper limit 
is taken as as a systematic error due to the line shape 
uncertainty of the J/'0-related background. This error 
is found to be 7.3% and 6 . 8 % for the XcO and Xc 2 se¬ 
lection conditions, respectively. By varying the order of 
the Chebyshev polynomial from the third to second, we 
obtained a set of upper limits from which we take the 
relative difference as systematic uncertainty due to the 
Xc0,2 —t TT^KgRir background line shape. The uncer¬ 
tainty is found to be 0.1% and 5.0% for the XcO and Xc 2 
selection conditions, respectively. By changing the mean 
of the Voigtian within 1 MeV/c^ and the width within 
1 MeV, thereby taking conservatively into account the 
uncertainty in the published mass and width of the rj^. [5] , 
sets of upper limits are obtained from which we take the 
maximum difference with the nominal upper limit as sys¬ 
tematic uncertainty due to the signal line shape. This 
uncertainty is found to be 2.3% and 2.5% for the XcO 
and Xc 2 selection conditions, respectively. The various 
systematic uncertainties on the fitting range and the line 
shape for the signal and background are highly corre¬ 
lated due to double counting of possible uncertainty con¬ 
tributions. The total systematic uncertainty on fitting 
is estimated by adding the individual systematic uncer¬ 
tainties in quadrature. The total fitting uncertainty is 
9.6% and 17.4% for the XcO and Xc 2 selection conditions, 
respectively. 

The systematic uncertainty of the number of ■0(3686) 
events is estimated to be 0.8% as reported in [17]. The 
uncertainty originating from the trigger efficiency is esti¬ 
mated to be 0.15% P5] . 

All the uncertainties on the branching fractions of the 
decaying particles of the channels of interest are obtained 
from the PDG |2] and are taken into account in the sys¬ 
tematic errors of our measurements. The corresponding 
values can be found in Tab. jTj 

Assuming that all the sources are independent, the to- 







tal systematic uncertainties 5 o ,2 are obtained by adding 
the individual uncertainties in quadrature. 

TABLE I. Summary of all considered systematic uncertainties 
(%). All uncertainties quoted are estimated to be symmetric. 


Source 

XcO Tr°r)c 

Xc2 ->■ 7r°77= 

Tracking of A, tt 

4.0 

4.0 

Kg reconstruction 

4.0 

4.0 

Photon reconstruction 

3.0 

3.0 

Kinematic 4C fitting 

1.2 

0.8 

PHSP generator model 

1.6 

2.7 

Four charged-particle tracks 

1.0 

1.0 

Ej resolution 

0.6 

0.6 

Fitting 

9.6 

17.4 

Number of 0(3686) 

0.8 

0.8 

Trigger 

0.2 

0.2 

B(0(3686) ^ 7Xcj) m 

2.7 

3.4 

B(77c ^ K^sKtv) id 

6.8 

6.8 

B(Ag ^ 7r+7r“) [2] 

0.1 

0.1 

Total So ,2 

13.8 

20.2 


V. RESULTS AND DISCUSSION 

The upper limits on the branching fractions of the 
XcJ 7 r° 77 c (J = 0, 2) transitions are calculated using: 

0 

B{XcJ -t TT 77c) < 'yXcj)Bintil - Sj)' 

where Nj^ are the upper limits on the number of sig¬ 
nal events, Sj is the total systematic uncertainty for the 
channel with J = 0,2, ej is the detection efficiency, 
Bint = B{r]^ ■ B{Kl tt+tt") • B(7r0 ^ 

77 ) = (1.7 ± 0.3) • 10“^ [2], and is the number of 

'0(3686) events [T7]. Table 0 summarizes the final re¬ 
sults of the analysis. 

TABLE 11. Summary of the final results for the XcJ (J = 0,2) 
decays. 



XcO TT°ric 

Xc2 ->■ 

vy® 

14.1 

35.9 

ej 

5.8% 

8.6% 

5j 

13.8% 

20.2% 

B{Xcj 

->• 7r°T7c)(10“®) < 1.6 

< 3.2 


In this paper, we presented an analysis with the aim 
to search for the hadronic isospin-violating transitions 
Xc 0,2 —t 7 r° 77 c using 106 x 10® 0(3686) events collected 


by BESIII through rjc —t decays. No statis¬ 

tically significant signal is observed and upper limits on 
the branching fractions for the processes Xc 0,2 
have been obtained. The results are B{xcO —t TT^r]c) < 
1.6 X 10“® and B{xc 2 —t < 3.2 x 10“®. These 

are the first upper limits that have been reported so 
far. These limits might help to constrain nonrelativistic 
field theories and provide insight in the role of charmed- 
meson loops to the various transitions in charmonium 
and charmonium-like states. Further developments in 
these theories will be necessary to clarify this aspect. 

The obtained upper limit on B{xcO does 

not contradict the theoretical estimate reported by 
Voloshin [5S] of order (few) x 10“'*. In this estimate, 
the branching fraction has been derived from a leading- 
order QCD expansion and related to the partial width 
of the decay 0(3686) —>■ h^TT^ under the assumption that 
the overlap integrals for the 2S—^IP and IP—J-IS tran¬ 
sitions are of similar value. In addition, Voloshin [29] 
predicts that the branching fractions of the hadronic de¬ 
cays Xco TT^rjc and Xci —t ' 7 r+ 7 r“ 77 c are approximately 
equal. A comparison of our result with that of an upper 
limit measurement of B(xci —t 7 r+ 7 r“? 7 c < 3.2 x 10“®) by 
BESIII [30] does not contradict such a prediction. We 
note, however, that an earlier theoretical estimate in the 
framework of a QCD multipole expansion m reported 
a branching fraction for Xci —t of (2.22±1.24)%, 

which contradicts the earlier BESIII measurement [30] 
and, under the assumption made by Voloshin |53|, our 
result as well. 

The near-future PANDA experiment |52| at the FAIR 
facility has the potential to find evidence or provide 
tighter constraints for the isospin-forbidden transitions 
discussed in this paper by directly populating the Xc 0,2 
states using an intense antiproton beam on a proton tar¬ 
get. 
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